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We numerically study the dynamical properties of fully 
frustrated models in 2 and 3 dimensions. The results ob- 
tained support the hypothesis that the percolation transition 
of the Kasteleyn-Fortuin clusters corresponds to the onset 
of stretched exponential autocorrelation functions in systems 
without disorder. This dynamical behavior may be due to the 
"large scale" effects of frustration, present below the perco- 
lation threshold. Moreover these results are consistent with 
the picture suggested by Campbell et al. in space of configu- 
rations. 
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I. INTRODUCTION 

At low temperature spin glasses (SGs) undergo a tran- 
sition characterized by the divergence of the nonlinear 
susceptibility. Moreover the relaxation functions of the 
system become non exponential already at temperatures 
higher than the transition temperature T sg . 

This behavior has been observed in canonical metallic 
and insulating spin glasses, investigated by neutron and 
hyperfine techniques 

In the Ising SG model, studied with spin-flip Monte 
Carlo dynamics, both in 2 dimensions (2d) ||] and 3d ||, 
non exponential relaxation functions have been observed 
below some temperature T* higher than T sg . Moreover 
in the 3d system Ogielski || observed that the long time 
regime of the relaxation functions is well approximated 
by the following function 



f(t) = f f 1 exp[-(t/ry 



(1) 



Fitting the data with this function he obtained that the 
onset of non exponential relaxation is consistent with the 
Griffiths temperature Tq, that coincides with the critical 
temperature of the ferromagnetic model. 

This result supports the argument suggested by Ran- 
deria et al. Q , recently verified by more rigorous analysis 
in Ref. || and by numerical simulations on a generaliza- 
tion of the SG model in Ref. (T^- According to these 
analysis in the SG the onset T* of the non exponential 
relaxation should be greater than or equal to the Grif- 
fiths temperature Tq. This behavior is caused by the 
existence of unfrustrated ferromagnetic-type clusters of 
interactions, the same that are responsible for the Grif- 
fiths sing ularity [§. The presence of non exponential 



relaxation in this approach is therefore a direct conse- 
quence of the quenched disorder. 

Another mechanism leading to non exponential relax- 
ation in frustrated systems, such as SG, has been sug- 
gested by several authors [p|-pd| . According to these ar- 
guments the onset T* of non exponential relaxation is 
greater than or equal to the percolation transition T p of 
the Kasteleyn-Fortuin and Coniglio-Klein (KFCK) clus- 
ters However, in frustrated systems with disorder 
T p is less than but close to Tq, therefore its eventual ef- 
fects are hidden by those related to Tq- 

A way to verify if percolation mechanisms can play 
a role in the dynamical transition of frustrated systems 
is to consider frustrated models without disorder where 
the Griffiths phase is not defined. In particular, we have 
considered fully frustrated (FF) spin systems |14| where 
ferromagnetic and antiferromagnetic interactions are dis- 
tributed in a regular way on the lattice, in such a way 
that no unfrustrated region (no Griffiths phase) exists, 
but the percolation temperature of KFCK clusters is still 
defined. 

In a previous paper B we have studied the 2d FF 
Ising model. We found numerically that the model ex- 
hibits a non exponential relaxation below the percola- 
tion temperature T p of the KFCK clusters. Moreover the 
long time regime of these functions is well approximated 
by a Kohlrausch- Williams- Watts function, also known as 
"stretched exponential" , 



f(t) = f exp[-(t/TY 



(2) 



In this paper we analyse, with conventional spin flip, 
the dynamical behavior of FF Ising model in 3d and in 
2c? using better statistics finding that T* is numerically 
consistent with T p . 

To clarify the role of percolation we study also the 3d 
g-bond FF percolation model. For q = 2 this model is 
obtained applying the KFCK cluster formalism to the 
FF Ising model (see Sect. ||). We simulate it using the 
"bond flip" dynamics 0. In this way the percolation 
properties of the model are stressed and the appearing 
of non exponential relaxation functions at T p are more 
evident. 

In both these cases we find that the relaxation func- 
tions exhibit an exponential long time behavior at high 
temperatures. Below the percolation temperature T p of 
the KFCK clusters, that is higher than the transition 
temperature T c of the model, the long time regime of the 
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relaxation functions becomes non exponential and is well 
approximated by a stretched exponential. 

Our results are consistent with the picture suggested 
by Campbell et al. || in the space of configurations and 
can be interpreted considering that T p corresponds to 
a thermodynamic transition in a generalized frustrated 
model § . 

In Sect, [n] we present the "g-bond frustrated percola- 
tion" model, and in Sect. [II we study the percolation 
properties of this model on a FF cubic lattice. We find 
that the percolation transition is in the same universality 
class of the g/2-state ferromagnetic Potts model confirm- 
ing the results obtained in the disordered version of the 
model in 2 d M . 

In Sect. |IV| we study the FF Ising model dynamical 
properties with conventional spin flip, and in Sect. fv| we 
present the relaxation functions obtained simulating the 
FF q-bond percolation model for q = 2, with the "bond 
flip" dynamics. 

In Sect.Jv3 we show th e co nnection with the Campbell 
scenario Qand in Sect. VIl we give the conclusions. 
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»'£[(• 

(ij) 



ijSiSj + i)S cri( j j 2], 



(5) 



II. THE "q-BOND FRUSTRATED 
PERCOLATION" MODEL 



in which every site carries two types of spin, namely an 
Ising spin and a Potts spin <Ji = 1, . . . , s with s = q/2. 
For q = 1 the factor q 1 *^ ) disappears from Eq. 
and we obtain a simpler model in which the bonds are 
randomly distributed under the conditions that the bond 
configurations do not contain a frustrated loop. For q — » 
we recover the tree percolation, in which all loops are 
forbidden, be they frustrated or not fT^ ]. 

When all the interactions are positive (i.e. = 1) the 
sum in Eq. (|4|) contains all bond configurations without 
any restriction. In this case the partition function co- 
incides with the partition function of the ferromagnetic 
q-state Potts model, which in the limit q = 1 gives the 
random bond percolation |l7j| . 

From renormalization group [ jl8| , mean field |l9[ | and 
numerical results f|@] we ex P ec t that the model in Eq. 
(||) exhibits two critical points: the first at a tempera- 
ture Tp(q), corresponding to the percolation of the bonds 
on the lattice, in the same universality class of the fer- 
romagnetic q/2-state Potts model; the other at a lower 
temperature T c {q) 1 in the same universality class of the 
FF Ising model. 



The FF Ising spin model is defined by the Hamiltonian 



n 



-'£< 



Eij Si Sj 



(3) 



where are quenched variables which assume the val- 
ues ±1. The ferromagnetic and antiferromagnetic inter- 
acti ons a re distributed in a regular way on the lattice (see 
Fig. [Vl| ). 

Using the KFCK cluster formalism for frustrated spin 
Hamiltonians |H| , it is possible to show that the partition 
function of the model Hamiltonian in Eq. (|^) is given by 



C 



(4) 



where q — 2 is the multiplicity of the spins, ks is the 
Boltzmann constant, \i = fcsTln(e 9j / fc£iT — 1), n(C) and 
N(C) are respectively the number of bonds and the num- 
ber of clusters in the bond configuration C . The summa- 
tion ~YT C extends over all the bond configurations that 
do not contain a "frustrated loop" , that is a closed path 
of bonds which contains an odd number of antiferromag- 
netic interactions. Note that there is only one parameter 
in the model, namely the temperature T, ranging from 
to oo. The parameter fi, that can assume positive or 
negative values, plays the role of a chemical potential. 

Varying q we obtain an entire class of models differing 
by the "multiplicity" of the spins, which we call the q- 
bond FF percolation model. More precisely, for a general 
value of q, the model can be obtained from a Hamiltonian 

M 



III. STATIC PROPERTIES 

In this Section we analyse the percolation properties 
of the model defined by Eq. (||) for q = 2, on a FF cubic 
lattice. 

After preliminary runs with spin flip dynamics on sys- 
tems with lattice sizes L = 10, 20, and with statistics 
of 5 x 10 3 thermalization steps and 5 x 10 6 acquisition 
steps, we found that the percolation transition occurs 
well above the critical temperature T c — 1.35 |21| (in the 
following the temperatures will be given in Jfk~B units). 
Then, we have simulated the model for L = 30 -f- 80, 
by the Swendsen and Wang cluster dynamics p^| , that 
turns out to be very efficient for the temperature regime 
of interest, allowing to consider only 5 x 10 4 acquisition 
steps. 

At every step we have evaluated the percolation prob- 
ability 



sn s 



and the mean cluster size 



s 2 n Sl 



(6) 



(7) 



where n s is the density of finite clusters of size s. 

Around the percolation temperature, the averaged 
quantities P(T) and S(T), for different values of the lat- 
tice size L, should obey to the finite size scaling p3] 
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P(T) = L- f) / v F P [L 1 / l '{T - T p )], 
S(T) = L->/ v F s [L 1 ' w {T - T p )], 



where /?, 7 and v are critical exponents, Fpix) and Fs(x) 
are universal functions of an adimcnsional quantity x. 

Standard scaling analysis results are summarized in 
Fig. |vn|. We obtained T p = 3.817 ± 0.005, v = 0.88 ± 
0.06, Pjt = 0.46±0.04 and 7/1/ = 2.03±0.03. The values 
of the critical exponents coincides, within the errors, with 
the random bond percolation exponents 

As we expect, the q =2-bond frustrated percolation 
model is in the same universality class of the q/2 =1- 
state ferromagnetic Potts model. 



IV. THE RELAXATION FUNCTIONS OF THE 
FULLY FRUSTRATED ISING SPIN MODEL 



(8a) random generator seeds, and between 10 5 and 10 6 steps 
for acquisition were taken, after about 10 steps for ther- 
malization. 

we show the relaxation functions obtained 
2.0, 1.8, 1.5, 1.0. For all temperatures we 



(8b) 



In Fig. 



VII 



for T = 2.5, 

fit the long time tail of the relaxation functions with Eq. 
(§■ 

The temperature dependence of exponents (3(T) is 
shown in Fig. VII. Note that (3(T) increases as func- 



tion of T from the value (3 = 0.61 ± 0.05 for T = 0.8 to 
the value = 1 for T > 2.0. 



As we can see in Fig. VII, our estimate of the onset 
of the stretched exponential relaxation is also consistent, 
within the errors, with the percolation temperature T p = 
1.701§. 

Within the errors the exponent x(T) increases as func- 
tion of T from the value x = 0.4 ± 0.2 for T = 0.8 to the 
value x = 1.6 ± 0.4 for T = 2.5. 



In this Section we present our results in the study of 
the FF Ising model, defined by the Hamiltonian in Eq. 
(|s|) for q = 2, simulated by spin flip dynamics. 

For each temperature T, 16 different runs were made, 
varying the random number generator seed, on a FF cu- 
bic lattice of size L = 30. We took about 10 4 steps for 
thermalization, and about 10 5 steps for acquisition, cal- 
culating at each step the energy E(t). The relaxation 
function of the energy is defined as 



fit) 



(SE(t)8E(0)} 



(9) 



where SE(t) = E(t) — (E). For each value of T, we 
averaged the 16 functions calculated and evaluated the 
error as a standard deviation of the mean. Here a unit 
of time is considered to be one Monte Carlo step, that is 
L d single spin update trials. 

In Fig. we show the results for T = 4.0, 3.5, 3.0, 



2.0, 1.5. We observe a two step decay also for high tem- 
peratures. For all the temperatures we fit the long time 
tail of the relaxation functions with the empirical formula 
proposed by Ogielski in Eq. (|l|). 

Temperature dependence of exponents f3(T) is pre- 
sented in Fig. VII. Note that (3(T) increases as function 
of T from the value (3 = 0.58 ± 0.03 for T = 1.5 to the 
value (3 = 1 for T = 3.7, 4.0. We do not observe any 
regular behavior in the temperature dependence of ex- 
ponent x(T). We estimated the errors on parameters as 
the range where we obtain a good fit of the relaxation 
function. 



As we can see in Fig. VII, these results are consistent, 
within the errors, with the scenario in which the onset of 
the stretched exponential relaxation coincides with the 
percolation temperature T p — 3.817 ± 0.005 (see Sect. 

D- 

We also simulated the FF Ising model on a square lat- 
tice of size L = 60. We calculated the relaxation func- 
tions of the energy. Averages were made over 16 different 



V. THE RELAXATION FUNCTIONS OF THE 
"g-BOND FRUSTRATED PERCOLATION" 
MODEL 

In this Section we analyse the dynamical behavior of 
the model defined by Eq. (ji|) with q = 2, simulated by 
the bond flip dynamics |7]]. 

The dynamics is carried out in the following way: at 
each step we choose at random a particular edge on the 
lattice; calculate the probability P of changing its state, 
that is of creating a bond if the edge is empty, and of 
destroying the bond if the edge is occupied; and, finally, 
we change the state of the edge with probability P. 

For each temperature T, 16 different runs were made, 
varying the random number generator seed, on a FF cu- 
bic lattice of size L = 30. We took about 10 3 steps for 
thermalization, and between 10 4 and 10 5 steps for acqui- 
sition, calculating at each step the density of bonds p(t). 
The relaxation function of the density of bonds is defined 
as 



fit) 



(Sp(t)S P (0)) 
<(M 2 > 



(10) 



where 8 pit) = pit) — (p) . For each value of T, we averaged 
the 16 functions calculated and evaluated the error as 
a standard deviation of the mean. We consider a unit 
of time to consist of G(p)^ 1 single update trials, where 
Q = 3L 3 i s the number of edges on the lattice. 

In Fig. VII we show the results obtained for temper- 
atures T = 4.0, 3.5, 3.0, 2.5. For T = 4.0, 3.5 we fitted 
the calculated points with the function in Eq. ([!]). 

The value of (3 extracted from the fit is equal to one 
within the error, and the value of x is zero. Thus for 
these temperatures the relaxation is purely exponential. 

For T < 3.5 we observe a two step decay, and only 
the long time regime of the relaxation functions could be 
fitted by Eq. (0). The value of [3 extracted is less than 
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one, showing that stretched exponential rela xation has 
appeared for these temperatures. In Fig. VII the values 
of f3(T) as function of the ratio T/T p are shown, with 
errors estimation. The exponent x(T) becomes non zero 
only for T — 2.5, for this value of temperature we obtain 
x = 1.1 ±0.1 



As we can see in Fig. VII, our estimate of the onset of 
stretched exponential relaxation is consistent, within the 
errors, with the percolation temperature T p = 3.817 ± 
0.005 of the KFCK clusters. 



VI. CONNECTION WITH THE RANDOM WALK 
PICTURE 

In this Section we make a connection between our 
model and the random walk picture of Campbell et al. 
Q , which we will briefly illustrate. Consider an hyper- 
cube in a D-dimensional space. Each summit is occupied 
with a probability p. On such a dilute lattice, a random 
walker is allowed to diffuse, like the "ant" on a percolat- 
ing cluster in the de Gennes picture. The mean square 
displacement after a time t is given by 



r 2 {t) 



(E£i(si(*)-si(0)) 2 ) 
D 



(11) 



where D is the hypercube dimension, x is a D- 
dimensional vector of components 0, 1, that identify the 
hypercube 2 D summits, and x(t) indicates the "ant" po- 
sition at the time t. 

Campbell et al. suggest in the Ising SG model that ac- 
cessible region in the space of configurations, compact at 
high temperature, becomes ramified at a temperature T* , 
and that a complex space of configurations is responsible 
for the appearing of non exponential relaxation. They 
also suppose that this temperature T* is the percolation 
temperature of the KFCK clusters. The idea is that the 
diffusive ant mimics quite well the evolution in the space 
of configurations in the SG model. 

In the study of the random walk on a randomly oc- 
cupied hypercube they find that for p < p* the function 
r 2 (t) becomes non exponential and is well approximated 
by a stretched exponential. But it is not possible to as- 
sociate a value of temperature to this probability. 

To make the connection between the bond frustrated 
percolation formalism and the random walk picture we 
introduce the local bond density autocorrelation function 

/(*) 



/(*) = 



£ti(Wi(Q))-(M*)) 2 
E?=i<&*(*)>-(M*)> 2 



(12) 



where Q = dL d is the number of lattice edges, d is the 
lattice dimension and bi = if the i-th bond is missing 
and bi = 1 if the i-th bond is present. The variables 
bi(t) are the coordinates in the space of configurations of 
the q— bond FF percolation model, which evolve by bond 



flip dynamics. In analogy with the picture proposed by 
Campbell et al., these variables can be interpreted as the 
coordinates of a walk on the summits of a ^-dimensional 
hypercube. For a fixed value of the temperature the walk 
will be confined in the subspace with density of bonds 
corresponding to that temperature^. 

Due to frustration not all configurations are allowed, 
therefore the walk occurs on a dilute space. By changing 
the temperature the space of configurations where the 
walk is confined changes and therefore the density of al- 
lowed sites in such a region of the space of configurations 
also changes. This is realized in a artificial way in the pic- 
ture made by Campbell et al., occupying randomly the 
hypercube summits (that represent the accessible states 
in the space of configurations). By changing the temper- 
ature, one may reach therefore a percolation threshold 
in the space of configurations. This would correspond to 
the breaking of ergodicity. At higher temperature, how- 
ever, the space of configurations may become ramified 
and stretched exponentials start to appear. 

Eq. ( |l2|) can be related to the distance travelled by the 
random walk r(t), via the relation r 2 (t) = 2{{p) - f(t)). 

We have simulated, by bond flip dynamics, the q = 2- 
bond FF percolation model on a square lattice of size 
L = 60. We have calculated the temperature dependence 
of the autocorrelation function in Eq. ([l2|). We find an 
exponential relaxation at high temperatures, while for 
T < T* the long time behavior of relaxation functions 
becomes non exponential and is well approximated by a 
stretched exponential. 

In Tab. | we show the fit parameters. Our estimate 
of the onset of stretched exponential relaxation functions 
is also consistent, within the errors, with the percolation 
temperatu re T v = 1.701. 

In Fig. VII we show the functions f (t) f or tempera- 
tures T = 3.5, 2.0, 1.7, 1.3, and in Fig. [Vl| the temper- 
ature dependence of the exponents (3{T) as a function of 
T/T p . 

We conclude that it is possible to apply the picture 
proposed by Campbell et al. to the g-bond frustrated 
percolation model. Furthermore our results are consis- 
tent with the hypothesis that the onset of non exponen- 
tial relaxation function coincides with T p . Note that we 
cannot exclude numerically that stretched exponentials 
are present even at temperatures higher than the per- 
colation transition, with an extremely small amplitude. 
This is also consistent with Campbell picture where the 
crossover from compact to ramified structure in the space 
of configurations is not sharp. 



The walk in such subspace is not random, since each step 
has a weight, which is a function of the clusters number vari- 
ation. If we consider the q =l-bond frustrated percolation 
model, the walk would be random. 
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VII. CONCLUSIONS 



In fully frustrated models we have numerically found 
an exponential relaxation above the percolation temper- 
ature T p , while for T < T p the long time tail of the relax- 
ation functions can be fitted with a stretched exponential 
in both 2d Q] and 3d systems. 

These results suggest that the percolation transition 
may play a role in the dynamical transition of frustrated 
systems without disorder. In particular, the role can be 
understood considering the physical meaning of the per- 
colation transition in a generalization of the SG model 
(the g-state Potts SG Q). We suggest that the percola- 
tion transition marks the appearing of the "large scale" 
effects of frustration. Below the percolation temperature, 
because of the presence of a spanning cluster, bond loops 
of any dimension may be closed, and therefore global ef- 
fects of frustration are observed. 

Note that in a previous paper we also studied 
a model, the "locally frustrated bond percolation" , in 
which only frustrated loops whose length is equal to four 
were forbidden. The model exhibits the same critical 
properties as the random bond percolation, showing that 
this kind of frustration is "too local" to change the uni- 
versality class of transition. Similarly, the relaxation 
functions in the long time regime can always be fitted 
with an exponential function. 
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T 


P 


r' 3 


5.0 


0.99 ±0.02 


0.16 ±0.01 


3.5 


1.00 ±0.02 


0.19 ±0.01 


2.5 


0.98 ± 0.02 


0.24 ±0.01 


2.0 


0.98 ± 0.02 


0.29 ±0.01 


1.7 


0.94 ± 0.02 


0.44 ±0.01 


1.3 


0.75 ± 0.02 


1.1 ±0.2 


0.7 


0.51 ±0.02 


2.8 ±0.2 



TABLE I. Fit parameters for the autocorrelation functions 
f(t) calculated with bond flip dynamics for the q — 2-bond 
FF percolation model 
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FIG. 1. Distribution of interactions for the FF model. 
Straight lines and wavy lines correspond, respectively, to 
Eij = 1 and €ij = — 1. 
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FIG. 3. Relaxation functions f(t) of energy as function of 
time t for the d = 3 FF Ising model, with spin flip dynam- 
ics, lattice size L — 30, for temperatures (from left to right) 
T = 4.0, 3.5, 3.0, 2.0, 1.5. 
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FIG. 2. Finite size scaling of (a) P(T) and (b) S{T), for 
the q = 2 model, and for lattice sizes L — 30,40, 50, 60, 70, 
80. 
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FIG. 4. Stretching exponents (3(T) as function of T/T p , 
the ratio of temperature over percolation temperature, for 
the d — 3 FF Ising model, with spin flip dynamics, lattice 
size L — 30. 
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FIG. 5. Relaxation functions f(t) of energy as function of 
time t for the d = 2 FF Ising model, with spin flip dynam- 
ics, lattice size L = 60, for temperatures (from left to right) 
T = 2.5, 2.0, 1.8, 1.5, 1.0. 



FIG. 7. Relaxation functions f(t) of bond density as func- 
tion of time t, for the q — 2 FF bond percolation model, on 
a d — 3 lattice of size L = 30, for temperatures (from left to 
right) T = 4.0, 3.5, 3.0, 2.5. 
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FIG. 6. Stretching exponents (3{T) as function of T/T p , 
the ratio of temperature over percolation temperature, for 
the d — 2 FF Ising model, with spin flip dynamics, lattice 
size L — 60. 



FIG. 8. Stretching exponents (3(T) as function of T/T p , 
the ratio of temperature over percolation temperature, for 
the q — 2 FF bond percolation model, on a d — 3 lattice of 
size L = 30 
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FIG. 9. Autocorrelation functions f(t) as function of time 
t for the q — 2 FF bond percolation model, on a d — 2 lattice 
of size L — 60, for temperatures (from left to right) T = 3.5, 
2.0, 1.7, 1.3. 



II 1 t 

0.9 - 
0.8 - 

; I 

0.7 - 
0.6 - 
0.5 - I 

0.4 ' 

0.5 1 1.5 2 2.5 3 

T/Tp 

FIG. 10. Stretching exponents f3(T) as function of T/T p , 
the ratio of temperature over percolation temperature, for 
the q — 2 FF bond percolation model, on a d — 2 lattice of 
size L = 60 
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